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ABSTRACT 


An automated digital computer technigue of ccntrol 
systes design is presented. The emphasis is on 
compensator design but the method is applicable to the 
design cf any system with free parameters. Signal 
representation and system response are in the time 


domain. 


Tke inputs. reguired from the engineer are the 
systen ccnofiguration, the desired output response and 
the free farameters. A parameter Dinimization routine 
is ther used tc minimize a specific cost function and 
to set the free parameters. A graphical cutput of the 
desired response and actual system response is tken 


produced fcr ccmparison by the engineer. 
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The past two decades have been witness to an ever 
increasing use cf the digital computer. Engineering usage 
and design are pretably the most important justificaticn for 
the large wmemcry, high speed computers of tcday. Classes in 
ccemputer aprlicaticn to engineering problems have teéccne 
part of tke established curricula at post universities. 
Numerous papers have appeared on the application of the 
computer to éngineering problems [1] which until recently 
were unsclvakle or at best required Vong, tedious 


precedures wkich aqave only approximate soluticns. 


Contrcl system engineering has relied increasingly on 
ccmputer sinulaticn of large scale systems to verify that 
design specifications have been met and ZO nake 
mcdificaticns to a system even before producing a protctype. 
There are pregrams available to simulate virtually any 
electrical cireuit, or control system, elther in transfer 
furction form cr as a system of first order differential 
equations. Cthers draw the Bode, Nichols cr Nyquist flots 
cf ecpen and clcsed locp systems. Some programs help design 
ccgfensatcrs which use an iterative methcd to achieve the 
desired frequency response [2] [3]- Researchers continue to 


better adapt the computer to engineering usage. 


Cantalapiedra [4] has used an iterative method tc find 
the optizum reduced order sodel for large order systems. 
MacNamara [5] went even further and used an iterative method 
tc find the cptimum compensation for an aircraft autofilot. 
It wculd eprear that these technigues cculd ke extendec and 


applied to the direct simulation and design of ccntrol 
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system circuits in the time dcmain. 


The irtert of this thesis was to develcp a user oriented 
program which cculd simulate a wide variety of ccntrol 
Systems and determine the values of the gains, feles and 
zeres necessary to preduce a desired response. A convenient 
means of data card input was to be provided to specify (1) 
the contrcl system which was to be optimized and (2) the 
desired reepense. A locally available function ninimization 
rceutine (ECXELX) was to be used to optimize the sinvlated 
system's cutrut. 


fo simulate the system which is tc be oftinized, 
transfer furctions which are commonly encountered were 
redauced tc first crder linear differential equations. These 
eqtations were then programmed so. that the transfer furction 
bliccks cculd be connected in an arbitrary fashion Ly data 
card infut. Several common nonlinear transfer blocks were 
alsc provided. The program will simulate the system with 
the kncwon parameters and then allow all unknown or 
adjustable farameters to he fixed by the ccgputer 


optimizaticn routine to achieve the desired response. 
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ri. PRCGRAM DEVELOPMENT AND IMELEMENTATION 


7k oe Gr ae Te < Sep Sr ae 2 2 oe ee ee > 2 Te me ee ae ae A owe see 


A. GENEFAL 


Any program which is to be of maximum tkenifit tc the 
user must kave a simple means of data input and an cutput 
which is e€asy to interpret and apply te the froblem at hand. 
The input data should have a physical significance that does 
nct lose its relevance through the programmirg of nunrerous 
egtuaticns. The program should be a readily usable tocl and 
not a proktler in itself. The intent of this thesis is _ to 
‘present suck a program, Computer Automated Lesign of Systems 
(CADS), which is readily used for Sinulation with 
Oftimizaticn cf the output. Optimum in the sense that the 
output resrense cf a given system configuration is as near 


the desired respcnse as possible. 


Most ccntrol system design starts with a prcposed 
transfer function block schematic to achieve a desirec time 
dcmain response. CADS has the ccmmon transfer functions 
built into tke pregram and the input data can come directly 
frcm the frorosed system schematic. The transfer tlocks 
which are availatle for system Simulaticn are presented in 
Tartle I. Tkese blecks should be adequate, either seperately 
or in varicus cascade combinations, to represent smost 


ccentrol systems. 
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To use the program the engineer must Know the systen 
configuraticn in transfer block form and the desired second 
order cutput response. If other than a second order 
response is desired, this may be easily specified but 
reguires scre knowledge of how to input the desired 
response. CALS will take the transfer block system, ccnnect 
it, ccmpare the given system response tc the desired 
response and set any free parameters to achieve the closest 


match possible of the system output to the desired restcnse. 


Figure 1 is representative of the type of system the 
prcegram car simulate and optimize. In figure 1, the 
parameters of the numbered transfer function blecks are 
knewn or fixed by e€guipment limitations. MTransfer blocks X, 
Y, and Z ccntain variable parameters which must be selected 
by the prcgram to make the system cutput rerrcduce a desired 


output furcticn as closely as fossible. 


- =~ 





FIGURE 1 Typical System to be Cfrtimized 
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E. MAIN FRCCEAM 


The MFAIN program was developed to contrel the selection 
of the sutfurcticns used in the optimization frocess and to 
ccgepute tke cesired response cf the systems which were to be 
cetimized. A second order ster response is calculated fron 
equation (I) as the desired response and stcred in the array 
XCATA. 

G 
@ - ses 2useve “Ct (I) 
The MAIN is essentially a bookkeeping routine which controls 
the Frcecraa execution. The subfuncticn operaticrs it 


centrols are defined in the following secticns. 


C. ELANT 


The ccmacn simtlation routines LISA, DSL, ECAP, CS¥*P and 
INTEG were irvesticated in an attempt to adapt them tc _ the 
ceneral frotlem specified above. These prcgrams required 
either the irput cf the system's state variable equaticns, 
Lapiace trarsform eéquations or the node-transfer function 
pair for each system simulated. These programs are useful 
for simulaticn if all of a system's parameters have been 
specified. Eawever, each simulation is problem specific and 
these teéechnigues are not readily integrated with cther 
sukfuncticn fregrars. A sSimulaticn subproaqram whickt was 
Capatle cf accepting data card input tc simuiate a wide 
Variety cf centrcel systems and of working with other 


existing subrrogrags had to te developed. 
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Subfurcticn PLANT was developed to simulate the systems 
which were tc ke optimized. It provides the versatility 
necessary to simulate numercus system configurations and is 
capable cf werking with a winimizaticn recutine to cptinmize 
the variakle ,arameters. The program reads the transfer 
furction blcck connections from data cards. This data is 
then used to autcmatically set up the system equatiors in 
state variakle format. These state variatle equations are 
sclved as first crder, ordinary differential equations by 
the Runge-Kttta-Gill forth order method. The programming 
was done in Fortran IV and all calculations are in dcuble 
precision. The capability of connecting the transfer 
furction biccks in any configuration (feed forward, 
feedback,etc.) .1s possible with this program. Simuiation 
flexibility is provided by having each Elcck within the 


system capakle of accepting an external forcing functicn. 


1, Eleck Data Card 


Several possible ways of inputting the data 
necessary tc simulate arbitrary system configuraticns were 
investigated. A means of defining a system configuration 
directly frem a bicck diagram schematic was chosen as _ the 
mest preferakle method. USing this approach, each data card 
was designed to be directly related to a transfer furction 
of the syste#. This resulted in having direct access to the 
transfer function farameters for optimizaticn. 


Ic simulate the system, a data card is prepared for 
€ach transfer bleck in the schematic diagram. The data card 
ccrctains a field cf numbers which specify the block ntnmker, 
tke type cf transfer function ccntained within the Elock, 
tke input ncde and the output node to which the Elicck is 
cconected and the values of any parameters associated with 


tre transfer function. The general fcrmat of the data 
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cards used tc input the system configuraticn is shown telcw. 


{1 faa {20 140 }60 (Column Numter) 
ELKCCLD=EFEVV G Pp Z 
Where: | 


CC = Pesiticn number of the block 
D 

EE 
VV = Cutfput ncde number 


Type of Elock (number) 


Input ncede number 


G, EF, and Z are parameters cf the transfer functicn. 


2. Eleck Connections 


The rumker of transfer blocks in the system and the 
Gata cards associated with each of the blecxs are read upon 
initial entry into the simulaticn routine. The program then 
connects the transfer function blocks in the proper order by 
cc@paring the input node number of a block to the cutput 
ncde nuagbter cf every other block. Whenever these two 
numbers are equal, the blocks are knowr to ke connected and 
a flag is set equal to 1.0. The flags are then used to 
identify the input drives to each block. 


The input quantity to each block is called the 
CRIVE. The program was designed to ailow for multiple 
inputs to the system teing simulated. This flexibility was 
acktieved by mgaking the input to each block equal to tre sum 
of the outputs cof all blocks connected to the input node 
rlus any external forcing function (DRVIN) feeding the 
transfer tlcck. The input DRIVE to a Elock is determined as 


shcwn in €guation (II). 


17 





DRIVE(i) = >> THA(3) *FLAG(35,i) + DRVIN(i) (II) 


maa) is the output of block j. FLAG(j,1) is 1.0 if block 
j is connected to Elock i and zero if it is not ccnnected. 
DRVIN(i) is any external forcing functicn specified Ly the 
user whick drives tlock i. DRVINS must be inserted in 
sukroutine ELANT as Fortran IV statements. The standard 
prcegram has CHVIN(1) = 1.0 specified as a unit step infut to 
the system. Proklem III-E in the section Investigation of 
Prcgram Ferfcrmance demonstrates how multiple, time varying 


inputs are tc bre inserted in the fprogran. 


Tre transfer function blocks available for system 
Simulaticn are shown in Table [. These blccks were selected 
because cf tkeir ccmmon usage in the modeling process. They 
were aisc fcurnd to be adeguate either separately cr in 
Various cascade combinations to represent most ccntrol 


systems. 


fre transfer function equations fcr each tyre of 
bicck were written in state variable format and stored in an 
array nagwed THA. The program reads a number from a block's 
Gata card which specifies the type of transfer function 
asscciated witk the block. The system eguations are then 
seclved sequentially by selecting the type of eqvation 
assccitaed with bicck number one, solving fcr its output and 
then sequencing tc block numker two, etc. The integral 
eguations are solved by a modified RKL fcurth order nrethod 
in the sukfunction rcutines. Three integration routines 
were necessary tc store the irtermediate results obtained 
for those equations involving a double integraticn. The 
fceur sutfunctions, RKLDE2, RKLDE3, CCPLX and RKIDE4 are 
called Ly ELANT. EKLDE2 is used for the integraticn of a 


18 





fiype two blicck. RKLDE3 is used for integration cf a type 
three blcck and to store intermediate quantities for 
Surruncticn CCPLX. CCPLX calls RKLDE3 and RKLDE4 for 
integraticn cf a type four block. 


Nc Erevision has peen made for the input of initial 
ccnditicrse cr the integrators. fherefore, the integrations 
must start at time egual zero. The user must specify the 
step size tc bre used for integration (DT) and the  feroklen 
run time (TIF). To conserve computer time [CT should be made 
as large as fossible. DT = TF/1000 is suggested as 
appropriate in most cases. Equally important is tnat TF be 
as small as fessible. Use of the program has shecwn that 
letting TIF te greater than the transient response time of 
the systen is rarely justified. For preliminary analysis fF 
was kept tc cnly Slightly lcnger than the time of the first 
undershoct fcr second order decminant systems. Because of 
the complexity of subroutine PLANT a flcw chart of PLANT has 


keen includec as Appendix B. 


fhe parameters of the system which are to he 
oftimized ky the mininization routine EOXPLX must be 
specified in PLANI. The minimization routine returns the 
trial values cf the variable parameters to ELANT in an array 
lakeled 'C'*. Regular Fortran IV statements equate the 


Variable parameters to the values in this array. If a 


pele-zerc fair were to he optimized, the fcllowing 
statements weoeuld be inserted in PLANT: P(i) = C(1) 
and Z(i) = C(2). The location fcr the preceding statements 


is clearly indicated in the program listing for PLANT. 
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D. DESIBED RESPONSE ({XDATA) 


The criteria against which the system response will be 
ccmpared will vary according to the applicaticn cf the 
system. Since most design work on control systems is done 
cn the tasis of second order dominance, the progrags was 
written tc provice a second order step response with 


adjustable 6, W and gain as the basic criteria against 
n 


which the sinulated system will be compared. The desired 6, 
.. and gain are read in as input data. The step respense is 
then ccmputed in the main program by subfunction RKLITEC and 
stored in the array called XDATA. Any cther time domain 
Lesponse nay be specified by the user by remcving the second 
ordér step response egquaticns and replacing them with the 
eguaticns of the desired response. An example is fgrcvided 
by problem III-E ir Investigaticn of Program Performance. 
If the pregram is Eeing used for simulation only and nc data 
curve is desired, setting the input variatle LEAP = 1 will 
causé the pregram to bypass the computation cf the standard 


seccend crder data equation. 


E. COSI FUNCTION (PERFORMANCE INLEX, FI) 


The achieved system resfonse is ccmpared witk the 
desired resjcnse and the difference is the errcer. The 
FrOgram sé€arckes for the parameter settings which will 
Minimize this error. The cost functicn may be specified by 
the user tc weight the system outputs as desired in 


sukfuncticn FE. The default cost function cf the program is 





TF 


tke integral e€rrcr sguared. Jj = (Err)? dt An exampl 
- orle 


of a weighted cost function is presented in Section III 
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i. MINIMIZATION RCUTINE 


1. seneral 


The free parameters are optimized to reduce the cost 
function Ey the complex methcd of M. J. Box.{6] Box's 
ccnstrained crtimization method has been programmed at the 
Naval Fostgraduate School by R. BR. Hilleary as a sukrcutine 
called ECXELX. This subroutine will find the minimum cf an 
acrkitrary function (cost function) subject to arbitrary 
explicit constraints and for implicit constraints. Explicit 
ccnstraints are defined as upper and lower bcunds on the 
free parameters. Implicit constraints may be arbitrary 
functicns of the free parameters (e.g. PP, <1768) . 


Iwc function subprograms are used to evaiuate the 
cbjyective functicn and implicit constraints, FE and KE 
respectively. The method BOXPLX uses te search fcr the 
vaiues of tke free parameters which minimize the cost 


function is explained in the ccmputer program listing. 


2. Explicit/Implicit Ccnstraint 


—_ <—- —_— <-_ =— 


1m 
\ 
ip 
ico 


ECXELX searches a feasibility region (n-dimensicnal 
Space, where n = number of free parameters) defined Ey upper 
and lower kEcunds on the free parameters for a minimus cost 
function. The sgaller the region defined by these 
koundaries, the amore rapidly the program will converge to 
the optigum parameter settings. Good engineering judcement 


will ke necessary to keep the feasibility region as smail as 


Pave 





possible. Tkte boundaries of the search regicn are read fron 
Gata cards Ey the main program as the upper bound (XU) and 
the lower bcturd (XL) for each free parameter. 


Implicit constraints may be any arkitrary furcticn 
of the free parameter desired. If an implicit constraint 
Such as tne product of a poie-zero pair must be less than 
scme numker is to ke evaluated, it must be supplied by the 
user to tke subfunction KE. No implicit constraints were 


used in tktis thesis. 


The starting values cf the free parameters (XS) are 
read in Ey tke MAIN frem data cards. A good chceice of 
Starting values wili dramatically reduce the time required 
fer optinizaticn. A preliminary root locus or Bode plot 
methcd of estimating the best values of the free variabies 


shculd ke accomplished whenever possible. 
G. GRAFHICAL OUTEUT 


Two sukroutines were written to provide for the 
graphical output of the desired response and the best system 
response achieved ry the optimization process. The user may 
select, by data card input, either subroutine PPL& waich 
prcevides a high speed printer flot or subroutine PIC which 
prevides a calcomp graph. Every fifth integraticn point 
stered in the arrays xXDATA and THAOUT ius flotted. 
Subroutines EFIT and FIC call the subrcutines PLOTP and CRAW 
respectively. FELOIP and DRAW are standard flotting routines 
at the NES computer facility and are not a part ci the 
Sinulaticn frogram. Figure 2 diagrams the information flow 


and data input to the progran. 
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FIGURE 2 cCADS Erogram Flow Caart 
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The example problems presented below were used to aid in 
the develcrment cf the optimization pregram. The order of 
difficulty cf the probiems progresses frcm a Simple text 
bock single variable, single input system tc a nultivariable 
Operaticnal servo drive system which has multiple inputs and 
discrete level feedback. An example of how a schematic 
diagram representaticn of a system iS prepared for infut to 
the program is presented pricr tc considering the example 


problems. 
A. CATA INEUT, AN EXAMPLE 


The Ward-leonard drive system [7] shown schematically in 
Figure 3 (a) has two variable parameters. The gain cf the 
amplifier and the tachometer feedback are available to 
adjust tke system's response. To simulate the system a 
blcck diacraz representation of the system is drawn as shown 


in Figure 3 (c) using the transfer blocks frem Table l. 
A = 
= ee 
ae ; —— ® 4 “oT 
— t ) | 


(a) 
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FIGCRE 3 Ward - Leonard Speed Ccntrcl Systen 
Using Feedback. (A) Schematic Ciagran 
(B) Elock Diagram (C) Bicck Diagram 
Using CADS Blocks. 


Where G(1) = K/r P(1) = Ifr 
G(2) = Km/yRa 
Gielen P(3) = 0 
G(4) = -Kn 
G(5) = -Kt 
DRVIN(1) = Er 
DRVIN (3) = -T, 

The nodes cf the block diagram are then nugbered 
Bermuentially 1,2,.-.<.«,0. The kElecks between the ncd€és are 
also nugbered seguentially ae a as shcwn in 
mragure 3 (Cc). Data cards (N) are then prepared for each 


bleck which specify the block number, type cf bicck, the 
input ncde, the output node, and the parameters ccntained 


Within tke kLicck. In Figure 3 (c), for example, biock 1 is 





a type two transfer block ccnnected Ltetween nodes 1 and 2. 
The data carc input for this block would be 


BL KO 12=0102 K/, i 


The program reeds the data card input and connects the 
biccks Ey setting a FLAG = 1. whenever the input node 
number to a Elock is the Same as the output nede number of 
any other tlcck. The input to a transfer block is then 
determined tc be the sum of the outputs of all tklccks 
ccnnected tc the input node plus any external fercing 
function driving the input node. The program has a unit 
step specified for DRVIN (1). If this is the only infut to 
the system, ne action is necessary on the part of the user. 
If cther than a unit step input to the system is desired or 
if there are other external forcing functions such as 
CRVIN (3), tkey must te specified and placed within the kody 
cf the sukfurction PLANT as Fortran IV statements. Fer the 
example shcwn in Figure 3 (c), a card with the equation 

DRVIN(3) = £(T ) 


would have tc be inserted preceding the drive equaticns. An 
example cf hcw multiple, time varying drives are specified 


1s given in section III. E,. 


When cptinizing a system, scme of the input quantities 
Will be unkrcwn or variables. These variables must also be 
assigned witkin subroutine PLANT. To optimize the variables 
C. and oe cf Figure 3 the follcwing two statements would be 


inserted in ELANT: 


G (1) 
G (5) 


C (4) 
C (2) 


where C(1) and cC(2) are the variables which will be 


optimized by subroutine BOXPLX. 
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Es TACHCMETER FEEDBACK 


The first optimization proktlem attempted was one which 
had an exact solution that can be found by algékraic 
methods. Arn instrument servo [8] with unity feedkack and 
fcerward transfer function 

1000 
S(S+10) (III) 


was to be cciapensated with tachometer feedback as shcwn in 


G (S) 


Figure YW. The cnly specification fcr the system's 
performance cf this single variable, seccnd crder system was 
the simple requirement that the closed lccp roots have a 
peee= 0.7. 






$@eo 
S$ (S+10) 







K,S 





FIGURE 4% Tachometer Compensated Systen 


The system shewn in Figure & was redrawn as Figure 5 in 
crder nc achieve the tachometer feedback. The 


characteristic equation of the system shown in figure 5 is 
S2 + (10 + 103 ‘ai + 103 = 0 (IV) 


Tke required Ss = 0.0343 may be calculated from egvaticn 
ay 
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FIGURE 5 cCADS Block Diagram of Tachcreter Systen 


CADS was frogramazed to optimize the system with the 
variable, ne specified we be between the limits 


O.C1 < a < 1.0. The desired response was specified tc _ be 
the standcarc seccnd order step respense for 6 = Q.7, 
Wn = 100C. CADS determined the optimum value for ss to be 
SS = 0.0343. Figure 6 shows the system's step respense and 


the desired response are virtually superimpesed. 


nS 





THAOUT 


0.8 


0.6 


0.05 0.10 C225 0.20 b225 
Time (sec) 


FIGUFE 6 CADS Compensated Tachcmeter Feedback 


System Response 
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The type one system, in the preceding example, allowed 
derivative feedback from the forward path without requiring 
a block which was capable of differentiation. One may 
enccunter a type zero system or a system which cannct be 
cenfigured to previde derivative feedback from the fcrward 
path. The pessibility of providing a fseudo derivative 
feedback using a type three bleck was investigated fcr these 


Cases. 


A type three bicck with Z = 0 and G = P is shcwn in 
Ficure 7. If the pole is placed far out orn the real axis, 
this blcck afrroximates derivative feedback. 





FIGUKE 7 Type Three Block used for 


Pseudo Tachncmeter Feedback 


The effect cf using this pseudo tachometer feedrack on 
the ccmplex roots of the closed loop system is neglicikle. 
It dees acd @ real roct at P = -964. 


The system of Figure 5 was redrawn as shown in Ficure 8 


using the pseudc tachometer feedback. 








FIGURE 8 Pseudo Tachometer Feedback 
Tre optinizaticn prcegram calculated SS = 0.0352 fcr the 


atcve configuration. This gave a 46 = 0.715. 


The system response and the desired response are shown 
mee Ligure S. The two responses are again rearly 
superinpecsed. This method of providing derivative feedkack 
has the disadvantage of adding an integraticn step tc the 
Eroblerm sciuticn with the concomitant increase in prctlen 


eclution time. 





THAOUT 


4 


FIGURE 


9 


AO 


0.05 0.10 0.15 0.20 0.25 
Time (sec) 


Fseudo Tachometer Ccmpensated System's Response 
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C. CASCADE CCMEFENSATION 


Having jfreven the feasibility of the prcegram optinizing 
a single variable system where an exact solution was 
availakle, tke next ysroblem considered extending the freoblen 
sccpe tc a third order plant with two variables. The 
unstable plant that was tc be compensated is shcwn in 
Figure 10. 









pee ee 
Cju)(jwrtC2jurn 


FIGURE 10 Third Order Plant to be Compensated 


The slant was tec be stabilized using a single section 
cascade ccmpensatcr. The compensated plant was reguired to 


have - <Z, wWithcut ceducing the error coefficient. {&] The 


compensated system is shown in Figure 11. Tc keep the error 
coefficiert ccnstant, the compensator used was a sirple lag 
network (7,< 7, ). 
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FIGURE 11 Lag Compensated System 


A ccupensator which would meet the recuired 
specificaticrs was calculated by conventional methods as a 
check or tke fregram's performance. The Bode pilot cf the 
uncompensated and conventionally compensated system is shcwn 
in Figure 12. The values of +7, and +r, fer Figure 11 which 
would meet the design requirements were determined tic he 
ec, = 1C., 7, = 100. The compensated system's respcnse 
using these values for the compensator is shcewn in Figure 


ie . 
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1.2 


eer 


THAOUT 


FIGUBE 13 


ZO 4.0 6450 8.0 308 
Time (sec) 


Corventionally Ccmpensated System Respense 
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The coupensated system was redrawn as Shcewn in 
Figure 14 using the standard program transfer tlocks 
available for system Simulation. 





FIGURE 14 CADS Lag Compensated Systen 


stardard second crder response of oo =a 
W = 90.75 was chcesen as the desired response for the 
n 
Optimizaticn frogram to match. The limits placed cn the 


cptimizaticn frogram were .001 < e <a Pan Gs, (20 e< - << te 
Arkitrary vaiues tc begin optimization were specified as 
FE = .07 anc 2 = .1 in the logarithmic centers of the 


Co Co 
search azcnes. CAES determined 7, = 8.078 and T, = 67.47 


as the ogtints farameter settings. Figure 15 shows the 
desired respense and the program ccmpensated systen 


response. 


The respense cf the system compensated by CADS is more 
nearly the desired response than is the conventicnaliy 
ccupensated systez. One should remember that the specified 
response is fer a true second order system whereas, the 


ccotpensated system is forth order. Therefore, a perfect 
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match cf desired versus actual responses could nct he 


cktained. 
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FIGGRE 1& cCADS Compensated Syster Response 
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rc. CASCACTE LEAL COMFENSATICN 


The ccmplexity of the next preblem to ke solved Ey CADS 
was extendec to five free parameters. The plant shcwn in 
figure 1€ was to be used tec follicw a unit amplitude 
Sine-wave input of 200 rad/sec. The output amplitude sas to 
re almcst exactly the same as the input amplitude, ané the 
output cculd not lag the input by mcre than 109. Two 


secticns cf cascade compensation were to be used. [8] 


—- ee 


_— 





FIGURE 10 Plant for Lead Compensation 


The specified requirements may be interpreted as an open 


lccp gain > 15 dk éand a phase angle ~ 90° at W = 200 from 
n 
a Nichols x_lct. A cut and try soluticr on a éEcde plot 


shewed tkat a gain of 3 XK 106 and tso phase lead 
comrpensatcrs with a double zero at Z = 70 and a double pcle 
at P = 7€C will satisfy the ciosed loop magnitude and fhase 
reguirements. Figure 17 shows the system response and 
desired respcnse obtained for these values. The magnitude 
of the ccampensated system's respense is 93% of the desired 


response and legs ry 8.12°. 
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The ercklem was then run on the optimization prcgran 
with the kleck connections as shown in Figure 18. The free 
parameters were the poles and zeros of the compensators and 
the gain of the plant. fhe desired data CULVE of 
XDATA = f€1r(z00t) was generated from the standard second 


order step response by setting 46 = 0, WH = 200 and using 
n 


X(z) as tke desired response. The initial searck zone 


Jinits were specified as 600 < P <800, 60 < Z <= 80, 
x a 


Bend 106 < ¢(3) < 3.2 x 106. 





FIGURE 18 CACDS Block Diagram for Lead Compensaticn 


The optigizaticn program sclution went to the lower 
linits fcr Ecth poles and the upper limits for toth the 
z€ros and tke gain. The limits of the search zones were 


relaxed tc 500 < Pp. <600, 80 2 eG ania 
a zi 


3.2 X 106 < €(3) < 3.4 X¥ 106. The program again places the 


free variakties cn the limits of P. = 500, 2 . = 9C, and 
nf ul 


G(3) = 3.4 X¥ 1C®. The system and desired response for these 
values is shcewn in Figure 19. The system magnituce and 
Ehase are much closer to the desired cutput response than 


tke respcnse cf the conventionally designed ccmpensator. 
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The phase cifference 1s only 4.339. Continued relaxation of 
the kouncaries produced the compensated system shcwn in 
Figure 2C. Figure 21 is a plot of the system's response for 
these values. The system's response is improved in that it 
lags the irput Ey oniy 3.469 and the magnitude 1s 
essentially tre same as the input magnitude. 


sinQoot-4#") 





FIGUEE 20 CALS Compensated System Elock Diagran 


The opftigizaticn program was activated at T = Q 
although tke prcktlem specifications were for the steady 
state resrcnsé. The problem was rerun with the optimization 
precess started after the transient had died out. The same 
values fer the cptinunm compensator were Crtained. 
Apparently the sazajll initial transient did not effect the 
prcklem scluticn. 
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ve SERVC SYSTEM COPPENSATIGN 


The fcregcing examples of the Simulaticn - optimizaticn 
Frcegram were simple examples which could obviously be 
guickly esclved with the standard cut and try metheés. A 
mcre complex and challenging example cf pregram usace is 
presented ky the system Shown in Figure 22. 


The system is an operaticnal servo drive mechanism with 
multiple infuts and discrete level feedback. This bighly 
ncenlinear system's output was to follow the input as closely 
as possikie and in steady state (t2 75ms) there was tec he 
very little reise ripple. The free parameters of the systen 


are the fecles ane zeros of the compensatcr, ous and the 


Ecoles cf the ncise suppressor C . 
n 


Tc sisulete and optimize the system it was redrawn using 
the availabie program blocks as shown in Figure 23. fupeng 
Situlaticr it was found that the current limiter for a. was 


net needed Lecause et < 25 A and the limiter was removed 


tc decreasé program run time. The desired response ({XDATA) 
was writter as a set of three eguatiocns. These equations 
were then used to replace the second order step respense 
eguaticns in the standard program. DRVIN(1), (4) and _ (5) 
were alsc writter as a set of egquaticns and placed in 
Subroutine ELANT. The discrete level feedback to bicck two 
was achieved by making DRVIN(2) = INTGER(THA(12)). These 
changes tc tte mair program and FLANT were all that were 
neccessary tec simulate this system. The insplementation of 


these changes to the program is shcwn on pages 80 and 81. 
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The system cptimizaticn was initially Eroken intc two 
parts. This was te reduce the number of variables that were 
to ke optimized per run and to obtain near cptimum starting 
values fcr the free farameters. The above steps were taken 
in an effort to decrease the computer time required for a 


scluticn. Tke first run was tc cptimize the compensator, 
eof independent of the noise supfresscr. The seccend 


Cptimizaticn cun was to select values for the noise 


SUPELESSCL, oe The rationale behind this separaticn was 
that the two circuits perform different functicns and should 


therefcre tbe initially separable in their effects cn the 
System. The final optimization run was to Fe made with ail 
free parareters available to the program fcr optimization. 


Tke search zcne centered on the values found above. 


Figure z4 shows the response for the syste as 


Criginally ccmpenseted. The oftimizaticn trun to set the 


values LCE S resulted in Z = 43.5, P = 21.0, 
Z| = 47.5, EF = 5$2.C0. Figure 25 shows the sisulaticn of 
the system tuSing these values. The initial velocity 


overshoot has keen reduced and the average velocity after 
the trarsient aprears more equally distributed akove and 


kelow the desired velocity. 
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fhe initial cptimization run for setting the valves of 


tke noise suppresscr used the same cost function whick had 


been used for all previous oftimization runs 

J = f ext dt. This resulted in ¢g¢ = 0.2 (lcwer bound), and 
W = 16CO (upper tound). Figure 26 shcws the result of 
D 


using this cost function was to further decrease the 


cvershoot. Ecwever, it allowed larger switching cr noise 
transients as a resuit of weighting large transient errors 
mcre than the lesser noise jitter. Tc overcome this 
prctlem, tke ccst function was changed to J -f [Err] *t*dt. 
This was tc weight the steady state errcrs mscre heavily than 
tke transient errcrs. Using this cost function the values 


set by tke cptimization prcgram were fd = .Z and H = 1430. 
n 
The damping factcr, d , was again placed on the lower limit. 


The initial cvyershoot was still imprceved cver the original 
System's resrpcnse tut the switching transients were not 
improved cver the previous optimizaticn trial. The results 
of the sinuletion run using these values is shown in figure 
a. 
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A final optimization run was made with all six 
variables free. The performance index was changed tc the 
weighted cost function shown in eguation V. The weichting 
factcr was ccmputed to equally weight the start up transient 
and the steacy state responses in an effort to reduce the 


Switching transients. 


J =f {Errjdt 0.0 < T < 0.045 

J = J4.14*[Erri dt 0.045 < t < TF (V) 
The optinta ccmpensatcr values established by CADS were 
P. = 7 Z = 52.5, PS =O 2 lisg Z = 48.1, 6 = 0.22 and 
a = 126€. Figure 28 shows the fully compensated system's 
response. The initial overshoot has been reduced and the 


average velccity in steady state is clcser tc the desired 
value tkan the criginal system's response. However, the 
transient errcr cn shut off is still present. Paro ie) sede 
summarizes tke results of the optimizaticn of the servo 


systen. 


The CEU time required for the optimizaticn precess was 
net decreased Ly spelitting the problem intc two seéfarate 
parts. The time required fcor each optimization run cn the 
individual ccrrensators was the same as the time reguired tc 


optimize tke ccmplete system with six variatles. 
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This frerlem has presented an excellent example cf the 
necessity tc carefully select the cost function which will 
Geasure the systea's performance. Defining a perfcrmance 
index which will weight the nore objecticnable 
characteristics of a system's response so that they are 
eliginatec¢ or reduced beccmes difficult as the system's 
ccuplexity ircreases. The performance index, J =/(Err)2 dt, 
was not adequate in its treatment of the ncise and switching 
trarsients fcr the above prcblen. The cther perfcrmance 
indexes tsed were also marginal in their effects upon 
parameter oftigization. Although the ccest functicr was 
reduced to a sathematically ccrrect tinipnun, the systen's 
perfcrmarnce tas not the best that could ce achieved. The 
system's perfcrmarce waS only optimum due tc the defirition 
of the ccst ftrcticn. The system was simulated usinc the 
ccmréensatcr values determined from the last cptimizaticn run 
with the excerticn of 6 which was increased to f= 
This simulaticn was made based upon engineering judgment of 
the effects cf varying ¢ . Figure 29 shows that this change 
in the damping factcr produced a system response which was 


nearer tke desired response than any of the frevious rtns. 
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A. CISCUSSICN 


The okjective of this thesis was to investigate the 
feasibility cf developing a computer program which wceuld 
optimize a variety of control systems' respcnses in the tine 
dcmain. he fregram develcped during the investication 
preves that time dcmain optimizaticn cf contrcl system 
responses is net only feasible but desirakle due tc the 
readily interpretable results of the optimization frecess. 
CACS requires appreximately 200K bytes cf ccuputer core when 
the high speed printer plot is used for graphical display of 
the output ard 230K when outputting calcomp gLapas. These 
ccre reguirenents are a maximum and could te reduced Fy more 


careful, prefessional programming. 


The ccuoputer time required for CALS to arrive at a 
scluticn is déefendent upon 

(1) the crder of the system being simulated, 

(2) the area cf the search zone determined by the 

upper and lower bounds on the variable farameters 

and 

(3) the nearness of the starting guess to the 

cptinup parameter values. 
Every trial value cf a parameter selected by FOXPLX requires 
a complete sirulation of the system in crder to evaluate the 
system's response and compare it with the desired resfense. 
The program run time is therefore, T,., = the number of 
trials xX system simulaticn time. A high crder system may 
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reguire twenty seccnds of CPU time for Sipnulation. If 300 
trials are reguired to determine the cptimum farameter 


values, tke tctal CPU time weuld be 100 minutes. 


Example Lower Starting Upper Number CPU Time 


Problem Bound Guess Bound or Required 
ET XL XS XU Trials 
P =300 -34.3 = 55  7min03sec 
(pseudo) 
~O01 2 Ou Oar 229 24min53sec 
Cc ~O1 1 1.0 225 24min53sec 
400 450 500 
400 450 500 
d 90 95 100 2,000 68min 
90 95 100 


3.4x10© 3.5106  3.5x106 





10 20 : 30 
40 50 60 
e 400 500 600 230 4nr 
aS 2 oD 
eS © S00 1400 
Table III 


fable III presents a summary of the search zcnes and 
tines reguired to cktain a solution for scme cf the frecblems 
ccnsidered ir this thesis. The amount of CFU time required 
fer a soluticn, e¢specially fcr problem III-E, may seen 
Excessive. Ecwever, there are several consideraticns which 
Smeald Ee mace pricr tc arriving at this ccnclusion. 

(1) The equaticns of the system do not have to be 


written, programmed nor detugged if the system's 
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compenent transfer functions are known. 

(2) Time domain requirements do not have to be 

translated irto frequency domain specificaticns 

for system simulation and design. 

(3) <A systematic search is carried cut to oktair 

the cptinum farameter settirgs. This assures that 

with valid -Eounds on the variables an acceptakle 

scluticn will be obtained with the Pipe t 

Ooptimizaticn attempt. 
The time reguired tc perform the above steps in the éesign 
precess Ey ccnventional means may result if wany more hours 
of CPU tine than if the program CADS were used fren the 


beginning of the design process. 


Several means of reducing the computer time recuired 
fcr cptinizaticn were previously outlined in secticn II. 
The most Significant reduction is obtained Ey keepine the 
Nugker cf integrations regquired for system Simulaticn to a 
Mirimum. Elcck diagram reduction of the system shculd be 
acccmplisted whenever possible. The example of the dard 
Lecnard Grive system shcwn on page 26 can ke reduced tc the 


Simple system shown in Figure 30 by blcck diagram reduction. 





FIGURE 3G Elock Reduced Ward~Lecnard Drive 
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Section IIJI.E. showed the effect of using different 
cost functicns tc measure the quality cf the cptigzized 
response. Alttough the integral error squared criteria is 
often used tc judge a system's performance, the user shculd 
carefully ccrsider how to best define a cost functicn which 
will preperly penalize deviations of the system response 
from the desired response. A small error (err < 1) scuared 
keccmes even smaller. If all the errors are small, the IES 
is a valid ccst functicn but if the system response invcives 
larce and sg#zall errors a weighted cost functicn will have to 
be used. One grethcd cf arriving at a properly weighted cost 
function is te simulate the system using first estimates of 
tke variakle farameters and recording the sur of the eérrers 
over the different portions of the respense. A ratio cf the 
errers can tken be used tc arrive at froper weichting 


factors fcr each time section cf the respense. 


BOXELX will ccntinue the optimization frecess, wcrking 
in the seventh significant digit until it can no lcnger 
reduce tke ccst function. Often an acceptarle soluticn for 
the system farameters has been found lorg before the ccst 
function kas teen reduced to its minimum. A judicious use 
cf CADS gay ke made by evaluating the system response after 
ten to fifteer minutes of run time to see if an accertable 


ecluticn has Leen found. 


B. SONCIUSICNS 


Time domain cptinization using the CADS program is a 
simple, straight forward process which dces not require an 
in-derpth analysis cf the system being optimized. Eccncmic 
use of CPU times does dictate that intelligent starting 
values and bcunds ke placed on the variable parameters which 


are to be optinizead. 





The frogram 1S a readily usable tool for simulation 
without cptinsization. It is ccmpetitive with, if not 
Superior to, other common Simulation routines when 
Situlatirng typical control systems. This feature alcne is 


expected to Iring the program into commcn usage. 


Cc. RECCHMENTATIONS FOR FUTURE WORK 


(1) All integral calculations rerformed during this 
investigaticn were done in dcuble precisicn for increased 
accuracy cf the system respcnse. The possibility cf using 
Single precisicn calculations should be investigated tc 


decrease ccre reguirenents. 


{Z) The akility to begin the optimization process at some 
time greeter than zero shculd be provided. This will 
necessitate ncdificaticn of the data input cards and block 


e€guations sc that initial conditions can be entered. 


(2) The program presently requires that external fcrcing 
furcticnse (CHVINS) and the variables which are to be 
optimized Fe specified by Fortran IV statements flaced 
within tke kEcdy cf the pregram. A method cf reading these 
specificaticrs frcem data card input should be develcped sc 
that the uwuser will not have to "shuffle" cards in the 


prcgram deck. 


(4) The craphical output of CADS was all that was necessary 
for the investigaticns conducted in this thesis [Fut a 
prcevisicn fcr numerical output of selected responses shculd 


ke previded fcr detailed analysis. 


(5) The feasibility of reducing the numker of significant 
digits BCXPLX uses should be studied as a seans of recucing 
cptimizaticn tine. Also some criteria might ke develcread to 
rezect a system's response before TF is reached if it is 


deterfined tc te unacceptable. 
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(6) A method of automatically relaxing the boundaries on the 
variables being optimized when they go to their limits 


should be developed. 


(7) The standard cost function provided and all user 
developed cost functions should be normalized. This would 
permit a more direct and easier comparison of a system's 
"goodness" when several different integration step sizes or 


run times have been used in the optimization process. 
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APPENDIX A 


Block Data Card Format 


11 114 120 440 #460 
ELKCCD=EEVV G P Z 

CC = PCSITICN NUMBER OF THE BLOCK 
[ = TYEE OF BLOCK (NUMBER) 

EF = INECI NCDE NUMBER 

VV = OUTEDOT NODE NUMBER 

G = VAIUE CF GAIN 

— = VALUE CF THE POLE* 

Z = WAICE CF THE ZERO* 


* For klicck type 4, 8 is read into the P locaticn and 


Wn is read into the Z pasition. 

* Fer kElock type 5, SF is read intc the P EQSLECIOR 
and or is read into the Z position. 

* EFcr Elcck tyfe 6, a is read inte the P fosition 


and ©@ is read into the Z position. 


EXAHFELES 
G 
PLKC 11=0102 10. 
G 
BL KC 22=0304 1. Be 
G ie 
BIKC23=0405 1. i eee 
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BLK 1C4=1006 eae .20 «14.14 
G¢ 8 -8 
BLK118=1207 en: 20beu=S: 
G c O, 
BLK126=0708 eee oer 1 0 


*See Table I. 
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APPENDIX B 
Plant Flow Chart 
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